Analysis of current-voltage characteristics of two-dimensional superconductors: 
finite-size scaling behavior in the vicinity of the Kosterlitz Thouless transition 
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It has been suggested [Pierson et al., Phys. Rev. B 60, 1309 (1999); Ammirata et al., Physica 
C 313, 225 (1999)] that for two-dimensional (2D) superconductors there exists a phase transition 
with the dynamic critical exponent z ~ 5.6. We perform simulations for the 2D resistively-shunted- 
junction model and compare the results with the experimental data in Repaci et al. obtained for 
an uhrathin YBCO sample [Phys. Rev. B 54, R9674 (1996)]. We then use a different method 
of analyzing dynamic scaling than in Pierson et at, and conclude that both the simulations and 
the experiments are consistent with a conventional Kosterlitz Thouless (KT) transition in the ther- 
modynamic limit for which z = 2. For finite systems, however, we find both in simulations and 
experiments that the change in the current-voltage (I-V) characteristics caused by the finite size 
shows a scaling property with an exponent a ~ 1/6, seemingly suggesting a vanishing resistance 
at a temperature for which z = a~^ . It is pointed out that the dynamic critical exponent found 
in Pierson et al. corresponds to the exponent a~^ . It is emphasized that this scaling property 
does not represent any true phase transition since in reality the resistance vanishes only at zero 
temperature. Nevertheless, the observed scaling behavior associated with a ~ 1/6 appears to be a 
common and intriguing feature for the finite size caused change in the I-V characteristics around 
the KT transition. 

PACS numbers; 74.76.-w, 74.25. Fy 74.40. -l-k, 75.40. Gb 
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I. INTRODUCTION 

The transition from the superconducting to the resis- 
tive state is of the Kosterlitz Thouless (KT) type for an 
ideal two-dimensional (2D) superconductor in zero mag- 
netic field.El Below the KT transition the current-voltage 
{I-V) characteristics is nonlinear, V oc 1°^ (or equiva- 
lently, E oc J"' with the electric field E and the current 
density J), characterized by the exponent a. According 
to the conventional extension of the KT theory, a = 3 
precisely at the transitioruayd a becomes larger as the 
temperature is decreased.Ero The scaling argument by 
Fisher, Fisher, and Huse (FFH) in Ref. ra leads to the 
connection a = 1 + z at the KT transition from which 
the value z = 2 is deduced for the dynamic critical expo- 
nent 2.13 These predictions of z = 2 and a = 3 at the KT 
transition have been confirmed through numerical sim- 
ulations for the Jattice Coulomb gas with Monte Carlo 
(MC) dynamicSjcl the Langevip-type molecular dynam- 
ics of Coulomb gas particles,B and the 2D XY nindel 
with resistively-shunted-junction (RSJ) dynamics .□u'tHl 
The value a — 3 has also been found to be consistent 
with numerous measurements.Ej 

From the above perspective, the recent claim by Pier- 
son et al. in Ref. [ll| that z has a much larger value at the 
resistive transition is very intriguing and is in apparent 
contradiction to the established conventional view. More 
specifically, Pierson et al. in Ref. O have applied a dy- 
namic FFH scaling approach (we call it Pierson scaling 
in this work; see Sec. O) to existing experimental data 
and obtained z ~ 6 in many 2D systems including theo- 
retical models like the lattice Coulomb gas (Ref. 0), ^He 
films, and an ultrathin YBCO sample (Repaci et al. in 



Ref. |2|). We in our analysis focus on the Repaci et al. 
data for which a particularly good scaling collapse was 
obtained by Pierson et al. 

Our approach is, as directly as possible, to compare 
simulation data for the 2D RSJ model obtained in this 
work with the experimental data obtained for a ultrathin 
YBCO film by Repaci et al. in Ref. n3. (For brevity, we 
from now on call the former the RSJ data and the latter 
the Repaci data, respectively.) In case of the RSJ data 
we conclude that the established conventional view holds 
with z = 2 at the KT transition. We then compare the 
RSJ data with the Repaci data and show that the latter 
exhibit a similar evidence of z = 2 at the KT transition. 
Both the RSJ and the Repaci data show nonvanishing 
resistances below the KT transition. In the former case 
this is due to finite-size effects and in the latter it is due 
to some length scale which introduces an effective cutoff 
in the logarithmic vortex interaction, e.g., the finite size 
of the sample, the perpendicular penetration depth, or 
a residual weak perpendicular magnetic field (see Ref. w 
for details); We in the following use the term finite-size 
effects to cover all thesepossibilities. In the original anal- 
ysis of the data in Ref. |l2|the authors attributed the finite 
resistance below the KT transition to such finite-size ef- 
fects, which is in accord with the striking similarity we 
find to the RSJ data. 

We also find that the Pierson scaling can be made to 
work with z « 6 for the RSJ data which is closely the 
same z as found by Pierson et al. for the Repaci data. In 
addition, we find, by analyzing the resistive tails of both 
data sets through the use of a new scaling method, that 
the finite-size caused resistance tails show strikingly sim- 
ilar scaling behaviors with an exponent a sa 1/6. How- 



ever, our conclusion is that this scahng does not signal 
any real resistive transition, but reflects finite-size effects 
which will eventually disappear in the thermodynamic 
limit. 

Strachan et al. (Ref. n3) have, fi'om a thorough analy- 
sis of the Repaci data, shown that the Pierson method is 
too flexible to allow any firm conclusion and that the ap- 
parent vanishing of the resistance is really a consequence 
of the finite-voltage resolution threshold in the experi- 
ment. Somewhat related conclusions have been reached 
by Brown in Ref. O in the context of the vortex glass 
transition. This is also in accord with our findings. Nev- 
ertheless our conclusion is that the Pierson scaling result 
is, at least partially, a reflection of an actual scaling like 
behavior associated with the flnite-sizc caused resistance 
tail in the vicinity of the KT transition. However, since 
the resistance tail vanishes only at zero temperature for 
a finite system, we conclude that there exists no resis- 
tive transition at a finite temperature and hence no new 
transition. A possible underlying reason for the scaling 
behavior is still an enigma which remains to be resolved 
in the future. 

The present paper is organized as follows: In Sec. || the 
scaling methods used to analyze the data are presented, 
and in Sec. Ill we describe and analyze the RSJ data we 



obtain from the simulations. The Repaci's YBCO data is 
reanalyzed and compared with the RSJ data in Sec. flV 
Finally we summarize and discuss our findings in Sec. 



II. SCALING METHODS 

The dynamic scaling relation proposed by FFH (Ref. |g) 
for a d-dimensional superconductor in the vicinity of the 
transition is given bytJ 



E = Je 



'-x±(A''-'/T), 



(1) 



where E, J, ^, z, and T are the electric field, the applied 
current density, the correlation length, the dynamic crit- 
ical exponent, and the temperature, respectively, and x± 
denotes the scaling function above (+) and below (— ) the 
transition. It should be kept in mind that Eq. (m can 
be directly used only in a system whose linear size L is 
much larger than ^ (L 3> ^). 

Specializing to two dimensions, d = 2, Eq. (m) can be 
recast into the form 



E 
1 



r^x±(J?/T), 



(2) 



which is the starting point in Pierson et alS3 It should 
be noted that the above FFH dynamic scaling form (g) 
is valid only for a phase transition when ^ is finite both 
above and below the transition temperature. Accord- 
ingly, one should note that Eq. (0) is not compatible with 
a system undergoing a KT transition because the whole 
low-temperature phase is "quasi" critical with-||-|= oo 
and is characterized by a line of fixed points J^'t£l This 



means that only the scaling functioii-j:^;+ above the KT 
transition has any clear justification.EZI 

In the present paper we study in Sec. |l| the 2D RSJ 
model with a finite size L, which, instead of ^, serves as 
the large length scale cutoff at and below the KT tran- 
sition. At the KT transition the finite-size scaling form 
deduced from Eq. (||) by substituting ^ by L takes the 
form:! 



E 

7 



L-^f{JL), 



(3) 



where the scaling function f{x) satisfies /(O) =const. 
and f{x) oc x^ for large x, corresponding to the finite- 
size induced resistance R oc L~^ without external current 
and E oc J^"*"^ in the large-current limit, respectively. 
The conventional view certifies that the nonlinear I-V 
exponent a ^%J^ oc J", is equal to 3 precisely at the 
KT transitioro'cl fijom which z = 2 follows by the scaling 
relation a = z + iB 

Since the low-temperature phase in the KT scenario 
is described by a line of fixed points, each temperature 
T below the KT transition is characterized by its own 
scaling function hT^x): 



YjLgi (r)Y 



^ = hT[JLgL{T) 



(4) 



where R = \iinj^o{E/J) oc L ^^"^^ [compare with 
Eq. (|)], hxiO) = 1, h{x) oc x^(^) for large x, and gLiT) 
is a function of at most T and L such that a finite limit 
function goc{T) exists in the large-L limit. 

The temperature-dependent dynamic critical exponent 
z{T) can be obtained from a simple scaling argument 
(See Ref. 10 for more details): The characteristic time 
T, which is related with R hy R '^ 1/r, is inversely pro- 
portional to the rate of unbinding vortex pairs over the 
boundary, and is given by the escape-over-barrier rate 



1 



oc nL exp 



1 



^TCG 



InL 



(5) 



where n is the average number density of vortices and 
the escape barrier for a pair to unbind over the bound- 
ary is determined by the effective logarithmic vortex in- 
teraction (l/e)lnL. Here e and T^'~^ are the dielectric 
constant and the Coulomb gas temperature, respectively 
(see Ref. M, and both can be calculated in equilibriuirL. 
From the definition of z(T) {t ^ U') and Eq. (g) we findcl 



z{T) = 



eTCG 



(6) 



Furthermore T'-^'^ oc T/pq{T), where po{T) is the density 
of Cooper pairs which decreases with increasing T, result-, 
ing in T'-^'^ monotonically decreasing to zero as T ^ O.EJ 
On the other hand, e(T) > 1 and remains close to unity 
in most of the low-temperature phase. It then follows 
from Eq. (O) that z{T) increases monotonically as T is 



lowered, reaching z{T) = cxo at T = 0. One may note 
that since the KT transition is characterized by the con- 
dition 1/eT'^^ = 4 (Refs. | and ^ Eq. (|) resuhs in the 
conventional established value z = 2 at the KT transi- 
tion. 

Since E/JR oc x^^^^ for large x [see Eq. (Q)] and 
z{T) depends on T [see Eq. (^)], E/JR cannot be de- 
scribed by a T-independent scaling function. On the 
other hand, when the scaling variable x = JEg^iT) has 
small values, the possibility that E/JR is described by 
a T-independent scaling function with the dominant T- 
dependence absorbed into x, cannot a priori be excluded. 
We find from our simulations of the 2D RSJ model that 
this indeed appears to be the case and one of the main 
points in the present paper is the apparent existence of 
the scaling form 



E 

Ir 



(^JLgLiT)^ 



(7) 



for smaller values of JLgL{T) [compare with Eq. (^]. 
What can be concluded about the function gL(T) on 
general grounds? First we note that the resistance 
R oc n{T)L^^^'^' for a finite system size vanishes only 
at zero temperature but has a finite value at any nonzero 
temperature. This follows since z(T) < oo at any 
nonzero T and n{T) is of the thermal activation form, 

i.e., n{T) ex e-l'=°""'l/^'"^ (Ref. 0), which is finite at any 
nonzero temperature. A finite R by Eq. (0) implies that 
< gL (T) < cxD and our simulations of the 2D RSJ model 
are entirely consistent with this expectation. If, on the 
other hand, R would vanish at a finite T this would have 
implied a diverging gL{T) at this temperature. 

The finite-size scaling we are discussing here is on the 
face of it entirely different from the scaling form given by 
Eq. (0), which is correct only in the thermodynamic limit 
^/L -^ and constitutes the basis and theoretical under- 
pinning for the Pierson method in Ref. O. We again em- 
phasize that this thermodynamic scaling is only justified 
as the KT transition is approached from above and has 
no bearing on the nonzero resistance below the KT tran- 
sition caused by the finite size of the sample. However, 
from the practical standpoint of analyzing procedures the 
Pierson method has a connection to our finite-size scal- 
ing form E/JR = h{JLgL{T)): The former essentially 
amounts to assuming that ^ in Eq. (0) is proportional 
to i?~", where a is a T-independent constant, leading 
to the scaling form E/JR = h{JR^°') [we have in this 
consideration ignored an additional weak unimportant T 
dependence in Eq. (||)]. Consequently, if ^^(r) = A^i?^" 
where A^ is a constant which may depend on i, then the 
Pierson scaling in a practical sense is of the same form as 
our finite-size scaling, in spite of the totally incompatible 
theoretical ground. So from this point of view the ques- 
tion is to what extent the relation gL{T) — AlR~°' is 
fulfilled. If we allow a to be a function of both tempera- 
ture and size then the relation gL{T) — A^i?^"^*^^^ just 
serves as a definition of a^ (T) because both g^ (T) and 



R{T) are positive functions. Thus, in the present finite- 
size context and from a practical standpoint, the Pierson 
scaling amounts to the statement that aL{T) extracted 
from experimental data appears to be a T-independent 
constant to a good approximation. 

The scaling form Eq. (Q) imphes that R{T) oc 
[T(7l(T')]~^ where the proportionality factor in gen- 
eral depends on temperature. Using the definition of 
aL{T), gL{T) = AlR~"^'^'^\ we get the connection 
[TgL(T)]"^(^)^(^) = BL{T)gL{T), where Bl{T) in gen- 
eral may depend on both T and T, which means that 



IngLiT) = — 



In- 



L 



aL{T)z{T) Bl{T) 



InBLiT). (8) 



The point with this relation is that if the dominant T- 
dependence of R comes from the factor (T(7i(T))~^(^^ 
then the T-dependence of B[^{T) is unimportant. Eq. (g) 
then leads to the requirement z(T) < l/aL{T) where the 
extreme condition z{T) = l/ai,(T) for a T > 0, if it oc- 
curred, would correspond to a diverging ghiT) and hence 
to a vanishing resistance. This means that if aL(T) for 
a fixed size L to a good approximation was a constant 
a, then the implication would be that the resistance van- 
ishes at the temperature for which z(T) = 1/a. As dis- 
cussed above this does not happen within the finite size 
KT scenario. But it could happen that a^lT) is a con- 
stant, i.e., aL{T) — a, over a limited temperature region 
and if all the available data was within this region one 
could be tempted to conclude that such a transition did 
really occur at some smaller T outside this region. As 
will be shown in Sec. Ill we do find such a limited region 
for the 2D RSJ model. 



III. ANALYSIS OF RSJ DATA 

The 2D L X L XY model on a square lattice is often 
used for studies of the KT transition. It is defined in 
terms of the Hamiltonian 



H 



-Ej y^ cos( 



A), 



(9) 



where Ej is the Josephson coupling strength, the sum- 
mation is over nearest-neighbor pairs, 9i is the phase an- 
gle at site i satisfying the periodicity: 9i = 9i+Lx = 
Oi+Ly- Here we use the fluctuating twist boundary con- 
dition (FTBC) to allow the phase difference across the 
whole system to fluctuate, and include the twist variable 
A = {Ax, Ay) in the Hamiltonian (^, where r^ is the 
unit vector from site i to site j; A = corresponds to 
the standard periodic boundary condition. In particular, 
FTBC is very efficient when the system is driven by an 
external current since the periodicity of phase variables 
di is preserved^ 

In the RSJ dynamics the local current is conserved at 
every site, which, determines the equations of motion for 
phase variables :EI 



i 



''3 / . 

fc 



dk - Tjk ■ '^) + V]k], (10) 



where the primed summation is over the four nearest 
neighbors of j, Gij is the lattice Green function for 2D 
square lattice, time t has been normalized in units of 
h/2eric with the critical current ic and the shunt resis- 
tance r of a single junction, and rjjk is the thermal noise 
current in units of ic on the link from the site j to fc. The 
equations of motion for the twist variables in the pres- 
ence of an external current density J in the x direction 
read 



A^ = 



^ Y, sin(0. - e, - A,) + 7?A. - J, (11) 



L2 



(y% 



^'--kT. 



j^ 2^ sin(6i, - Oj - Ay) + rjAy, 



(12) 



which have been obtained from the condition that the 
net global current across the system in each direction 
should vanish (see Refs. || and |o| for discussions). Here 
^/j 1 denotes the summation over all bonds in the x 
direction, J is in units of ic (the lattice constant is set to 
unity). The thermal noise terms rjij in Eq. (IIQ), t/Aj,. in 
Eq. (11^) , and 77a„ in Eq. ( p^ ) obey the conditions {rnj ) = 

(vaJ^ ivAy) = inijVAj^ ivijVAy) = {■nA^VAy) = 0, 

{vaMvaA^)) - {VAyithAM) - (2T/L2)<5(i), and 
{ili3{t)'>lki{0)) = 2T5{t){5ik5ji-5ii5jk), where T is in units 
of Ej/ks- The voltage drop V across the system in the x 
direction is written as ^ = —LA^ in units of icr, and we 
measure time- averaged (denoted by (•••)*) electric field 
E = {V/L)t to obtain the I-V characteristics. The set 
of equations (pi)|), (|ll|), and ([l^ ) are integrated by using 
the second order Runge-Kutta-Helfand-Greenside algo- 
rithnxJ with discrete time step At = 0.05 for L = 4, 6, 
and 8, and the electric fields £"s are measured from time- 
averages over O(IO^) time steps for large currents and 
O(IO^) steps for small currents. 

Figure |l| shows the I-V characteristics (£' versus J) at 
various temperatures in log scales for the system size L — 
8. In the small-current regime E is shown to approach 
the Ohmic behavior {E/ J = const.) at any tempera- 
ture used in the simulations. Since it is well established 
that the 2D XY model on the square lattice has the KT 
transition at Tkt ~ 0.892 (Ref. |2^), the ubiquitous non- 
vanishing resistance tails in the small-current regime are 
clearly due to the finite-size effects, and in the thermody- 
namic limit the resistance R = \m\j^o{E/ J) vanishes in 
the whole low-temperature phase below Tkt giving rise 
to a nonlinear I-V characteristics: E ex J^(^)+i where 
z{T) = 1/eT*-^*^— 2|Can be readily obtained from the equi- 
librium propertiesJj At the KT transition z{Tkt) — 2 
which corresponds to £' oc J"^ (denoted by thick full line 
in Fig. 1^ at T = 0.90 very close to Tkt) and crosses over 
to the Ohmic behavior E/J =const. as J is decreased 
due to the finite-size effect (thick broken line). 



The nonlinear I-V characteristics for the somewhat 
larger currents is caused by the breaking of vortex- 
antivortex pairs due to the finite current. It is this pair 
breaking which is described by the exponent a(T) — 
z(T) -|- 1; it continuously increases from 1 in the high^ 
temperature limit to infinity at T = as T is decreased.cl 
Consequently one should note that there is no jump 
in the exponent describing the current induced vortex- 
antivortex pair breaking at the KT transition. We stress 
this point because this fact seems sometimes to be mis- 
understood and to cause confusions. 

Suppose now that the I-V characteristics in Fig. [^ was 
all information we had. What feature in the data then 
implies that we are looking at a KT transition broadened 
by finite-size effects? Figure |2| shows the same data as 
d In E/d In J plotted against J in log scale. Below the 
KT transition in the low-temperature phase the slope 
d In E/d In J for smaller J should be a constant given 
by z{T) -\- 1. However, because of the finite-size effect, 
the slope instead crosses over to 1 for smaller J. This 
crossover produces a maximum for the curves in Fig. g. 
Since the crossover is a finite-size effect this means that 
the position of the maximum is controlled by the variable 
JL up to additional T-dependencies. If the additional 
T-dependencies are weak this implies that the maximum 
is occurring at roughly the same /(= JL) for all T be- 
low the KT transition. This turns out to be true for 
the 2D RSJ model as indicated by the full drawn lines 
connecting the maxima in Fig. 0. Above the KT transi- 
tion the maximum is instead controlled by the variable 
JS, and since the correlation length ^ strongly decreases 
with increasing T above the KT transition this means 
that the maximum above the KT transition should move 
to larger J with increasing T. The crossover between 
these two behaviors should occur near the temperature 
for which L « ^, i.e., a T somewhat larger than Tkt- 
In Fig. H such a crossover occurs between T — 1.00 and 
1.10 which is consistent with Tkt ~ 0.90. The inset in 
Fig. y compares these maxima (denoted as 'max') with 
the values z{T) -f 1 expected in the thermodynamic limit: 
the one obtained from a static calculation [see Eq. (|g)] 
(denoted as 'static') and the other one from the scaling 
form R ~ L~^^'^'> (denoted as 'scale', data from Ref. |^) 
[see Eq. (0))]. As seen from the inset, the maxima in 
Fig. g agree very well with the actual values for the ther- 
modynamic limit. This implies that the KT transition 
occurs close to the temperature for which the maximum 
is 3, which according to the inset gives a value close to 
T — 0.90. This suggests that the signature of a KT tran- 
sition for a finite-size system is refiected in the behavior 
of the positions of the maxima and that the transition is 
close to the T for which the maximum is 3. 

Next we demonstrate how the usual finite-size scal- 
ing given by Eq. (^) works for the case of the 2D RSJ 
model. First we assume that the L dependence of ql (T) 
is weak, so that the scaling variable is to a good approxi- 
mation JL. FiguijeiSa) then verifies the expected result 
that z{Tkt) = 2sM This is done by calculating the /- 



V characteristics for the known KT transition tempera- 
ture Tkt ~ 0.90 of the 2D XY model for some lattice 
sizes and then plotting L{E/ J)^/^ against LJ. The z 
value is then determined by finding the best data col- 
lapse for small values of JL. Since the low-temperature 
phase is a line of fixed points, the same procedure can be 
used to determine z{T) for any T in the low-temperature 
phase. Figure y(b) gives the same scaling determination 
at T = 0.80 in which case the value z{T) = 3.3 which 
should be compared to the value z{T) = 3.46 obtained 
from the direct static calculation of z{T) — l/eT'~"~^ — 1. 
The determination in Fig. @(c) includes the L dependence 
of gL {T) [the determination of gL (T) will be given below] 
and with this included a good collapse is obtained for the 
value z{T) = 1/eT'^'^ - 1. Figure | thus illustrates that 
the exponent z{T) can be determined from the finite-size 
scaling, depends on temperature, and has the expected 
value z{T) = 1/eT^^ - 1. 

As discussed in Sec. ||, the existence of the standard 
finite-size scaling, as illustrated in Fig. 0, does not rule 
out the existence of a scaling of the form Eq. (|^) for small 
values of LJg]^{T). This scaling form implies that E/ JR 
is a function only of LJgL{T). In Fig. |4| E/JR is plot- 
ted against LJgL{T) and shows that a function gL{T), 
which makes the curves collapse, can indeed be found. 
The function gL (T) determined from the two sizes L = 6 
and 8 are shown in Fig. and in Fig. H gh^T) ior L — 8 
is plotted against z{T). We emphasize that the existence 
of a function g^lT), which makes the data collapse, is 
by no means trivial and that as a consequence the func- 
tions are well defined. We again stress that the scaling in 
Fig. ^ [with gL (T) in Fig. ^ is compatible with the usual 
finite-size scaling given by Eq. (H). 

As discussed in Sec. || the Pierson scaling analysis, 
which presumes a thermodynamic limit, is completely in- 
congruent with the finite-size scaling of the 2D XY model 
discussed here. Nevertheless a connection from a prag- 
matic standpoint would arise provided gL{T) = A^R^"', 
as also discussed in Sec. |^ Figure g demonstrates that 
such a relation does exist over a limited T region. The 
exponent a « 1/6 is found both for lattice sizes 6 and 
8 and is well determined since R has a very strong T 
dependence. As discussed in Sec. |^ in connection with 
Eq. (|), the relation gL{T) = AlR'°', if it applied aU 
the way down to the temperature for which z(T) = 6, 
would imply that R vanishes at this T. This is, from 
our finite-size perspective, the connection to the Pierson 
scaling method and the putative thermodynamic tran- 
sition. We illustrate this in Fig. 0(a) where the data 
for which the relation gL{T) = A^R^" holds is plotted 
in the Pierson form using the T^^ — 0.125 for which 
z{T) == 1/T'^'^ -2 = 6 which is where the would-be 
transition should occur. The connection is that 1/z in 
Fig. M(a) is our a and that the variable ^ corresponds to 
R~^'^ . Furthermore R^^^^ is assumed to have the form 
exp(6/^rCG _ TCG) where Tf'^' is where the would- 
be transition should occur. As seen from the inset, R 
can be represented by the suggested form with the "cor- 



rect" T^"^ = 0.125 and a good scaling fit is obtained for 
the value z = 1/a = 6. Again we stress that there is 
no real phase transition or a vanishing resistance at this 
T^'^; the resistance is caused by the finite size of the 
system and the actual scaling property is only present 
over a limited range. Indeed when the lower tempera- 
tures T — 0.70 and 0.65 are included, a Pierson scaling 
plot with a ~ 1/6 can no longer be obtained. How- 
ever, a Pierson scaling plot can again be obtained if a is 
changed to a « 1/12 and the two highest temperatures 
are excluded as shown in Fig. |^(b). This suggests .tiiat 
the Pierson scaling method has too much flcxibilityJl3 



IV. RE-ANALYSIS OF REPACI'S YBCO DATA 



In this section we compare our analysis of the 2D XY 
model with the measured I-V characteristics for an ul- 
trathin YBCO sample by Repaci et a/Jla Figure |8| shows 
Repaci's I-V data for lOK < T < 40K (from left to 
right) which should be compared to Fig. ^ for the 2D 
RSJ model. The thick lines have the slopes one, three, 
and seven (from left to right), respectively, where the 
slope one represents the linear resistive tail, slope three 
is the value for the KT transition and slope seven the 
value where the scaling over a limited region for the 2D 
RSJ model would place a nonexistent "ghost" transition. 
Repaci et al. attributed the linear resistance below the 
KT transition slope three with free vortices induced by 
the finite size of the sample (width 200 ^m). The sim- 
ilarity between Fig. y and Fig. g is in accord with this 
interpretation. Figure shows the Repaci data plotted in 
the same way as the data for the 2D RSJ model in Fig. ||. 
The similarity is again striking. As in Fig. the maxima 
in Fig. ^ occur at approximately the same / for lower 
T (maxima larger than 2). This is consistent with the 
finite-size interpretation given in connection with Fig. ^ 
that the position of the maxima is roughly determined 
by JL. Consequently, when the size L is the dominant 
length, the maxima occur at the same /. However, for 
higher T (maxima lower than 2 in Fig. S) the correla- 
tion length ^ is the dominant length (^ < L) and, since 
S, rapidly decreases with increasing T, the maxima are 
displaced towards higher /. In parallel with Fig. g the 
thermodynamic KT transition is expected to occur close 
to the maximum value 3 which means T w 27K (compare 
full drawn line with slope 3 in Fig. nf). Consequently ^ 
is expected to become the dominant length for a some- 
what higher T and this is consistent with the fact that 
the maxima seem to start to move roughly around max- 
imum value 2. The similarities between Figs. |l| and g for 
the 2D RSJ model and Figs, g and|^ for the Repaci data 
support the interpretation of the Repaci data in terms of 
a finite-size induced resistance below the thermodynamic 
KT transitions. 

The next question is then whether the size scaling 
Eq. (|^) for the 2D RSJ model likewise applies to Re- 



capi's YBCO data. Figure |o|(a) corresponds to Fig. |j(a) 
for the RSJ model. The data in Fig. |l^(a) represents 
the data in the interval for which a linear resistance 
V/ 1 = R IS clearly present in the measurements (which 
means 23K < T < 40K). The data collapse in Fig. |l^(a) 
shows that V/IR is only a function of the scaling variable 
JLgLiT) when JLgLiT) is small enough, while Fig. |l^(b) 
shows that the form of the obtained scaling function is 
very similar to the one found for the 2D RSJ model. This 
suggests that the cause of the scaling is the same for the 
two cases, supporting the finite-size interpretation. The 
obtained scaling function QhiT) is plotted in Fig. |ll|, and, 
just as in the 2D RSJ case, it is very well represented by 
the AR~°' with a = 1/6. The range over which this is 
true for the 2D RSJ model is roughly 2 < ziT) < 3.5 
(see Fig. ||) and the gL{T) in Fig. ^ in fact covers the 
somewhat larger range 2 < z{T) < 4.5. As explained in 
Sec. [I], it is the fact that gL{T) is well represented by 
AR^"' with a = 1/6 over a range around the KT tran- 
sition which singles out the temperature where z(T) « 6 
as the "ghost" transition. The temperature correspond- 
ing to the slope z{T) + 1 = 7 is given by the steepest full 
drawn line in Fig. pi However, in analogy with the 2D 
RSJ model, we infer that there is in reality no vanishing 
resistance except at T = 0. In this sense it represents a 
"ghost" transition. 

Since gLiT), a, and z{T) (estimated from the max- 
ima in Fig. pi) are known, we can directly test Eq. (0) by 
plotting In g^ (T) as a function of a2:(T)/[l — az(r)]. Ac- 
cording to the discussion of Eq. ^ this should roughly 
be a straight line. As seen from Fig. |lj this is indeed 
the case for T below z « 2 (marked by a vertical dashed 
line in Fig. 12, note that z k, 2 is approximately where 
the finite-size effect sets in according to Figs. and 0). 
If it continued to be a straight line for lower T then the 
condition 1 — az(T) = implies g^ ^ oo and R = 
signaling the "ghost" transition for z{T) — 1/a. 

An alternative way of demonstrating the existence of 
the scaling Eq. (M) with the scaling variable IR~°' is to 
note that this scaling implies that V/IR for a given IR^" 
is independent of T (see Ref. 21, where a similar method 
has been used). Thus we can pick a given V/IR and 
mark the value / and RI for which it occurs. As seen 



in Fig. 12 this construction gives rise to a straight line 
for a given value of V/IR in the plot of V against / in 
log-log scale and the slope of such a line is according to 
the scaling assumption l/a + l. As apparent from Fig. |l^ 
this alternative test of the scaling is also borne out and 
gives values in the vicinity of 1/a k, 6. 

Finally in Fig. ^ we plot the result from the Pierson 
scaling method for comparison. 



V. DISCUSSION AND CONCLUSIONS 

We have shown from simulations that the finite-size 
induced features of the I-V characteristics for the finite- 



size 2D RSJ model obey a scaling of the form given by 
Eq. (|^) over a limited temperature range in the vicinity 
of the KT transition. The T-dependence of this scaling is 
absorbed in a function gL{T) for each fixed size L. Fur- 
thermore it was shown that this function gL{T) over a 
limited region in the vicinity of the KT transition is to 
a good approximation proportional to R~°' . For the two 
sizes we studied (L = 6 and 8) we found that ghiT) was 
different, yet the exponent a w 1/6 was found in both 
cases. We noted that if the scaling form Eq. (0) was valid 
down to low enough T and that at the same time the 
proportionality gL{T) oc i?~" was valid to low enough T, 
then this would imply that the resistance vanishes at the 
T for which z{T) = 1/a. However, we stressed that in 
practice there is no such transition for the finite-size 2D 
RSJ model and introduced the term "ghost" transition 
for the transition which is not there. The fact that the 
transition does not occur is linked to the fact that the 
scaling in practice breaks down for lower T, as clearly 
shown in Fig. H. For the finite-size 2D RSJ model the 
resistance vanishes only for T = 0. 

We also stressed that the theoretical basis for the Pier- 
son scaling method is incompatible with the finite-size 
effect we study for the 2D RSJ model. However, in 
practice the Pierson method assumes that for the high- 
temperature branch of the scaling form Eq. ^ one has 
f ex i?~". Thus the high-temperature branch of the Pier- 
son scaling is, from a pragmatic point of view, connected 
to our finite-size scaling. We explicitly demonstrated this 
connection in Fig. for the finite-size 2D RSJ model. 

We compared the data obtained for the finitCpSize 2D 
RSJ model with the YBCO data by Repaci et a/.ll3 From 
the strong similarity manifested by Figs. H, 0, H, and O we 
concluded that Recapi's data correspond to a KT tran- 
sition around 27 K and that the resistance below this T 
is finite-size induced. This Js in accord with the original 
conclusion by Repaci et alx/n Here we want to emphasize 
that finite size essentially means a cutoff in the logarith- 
mic vortex interaction. Apart from the actual finite size 
of the system this can be caused by a small remnant 
magnetic field or by a finite perpendicular penetration 
depth.tl 

We demonstrated the existence of a scaling of the form 
Eq. (1^) also for Repaci's YBCO data and furthermore 
demonstrated that the shape of the scaling function ap- 
peared to be the same as for the finite-size 2D RSJ model. 
This we interpreted as that the scaling for the YBCO 
data and for the 2D RSJ model have a common origin. 
This further supports the finite-size interpretation of the 
YBCO data. The function gL{T) was just as in the 2D 
RSJ model well represented by gL{T) oc i?~" with the 
same exponent a « 1/6 as for the 2D RSJ model. The re- 
markable similarity between the results for the finite-size 
2D RSJ model and Repaci's YBCO data clearly suggests 
that the same physics is reflected in the two systems. 
Our conclusion is then that, since it reflects a finite-size 
induced feature for the 2D RSJ model, it also reflects a 
finite-size induced feature for the YBCO data. Conse- 



quently, we conclude that the transition to a vanishing 
resistance imphed by the condition z(T) — 1/a = 6 is 
also for the YBCO data a "ghost" transition. 

Pierson et al. in Ref. O, when analyzing the Repaci's 
YBCO data with the Pierson scaling method, inferred a 
transition to a vanishing resistance for z{T) = 5.9. This 
is closely the same value as where our finite-size scaling 
places the "ghost" transition both for the same experi- 
mental data and for the 2D RS J model. This we believe is 
because the high-temperature branch of the Pierson scal- 
ing from a pragmatic point reflects the same behavior as 
our finite-size scaling. The low-temperature branch of the 
Pierson scaling is as seen in Fig. nj less convincing and it 
does not have any relation to our finite-size scaling. The 
fact that Pierson et al. for very many cases have found 
that their scaling suggests a transition at z{T) w 6 from 
our perspective suggest a common origin. The fact that 
the same result is consistent with the finite-size 2D RSJ 
model suggests to us that this common origin is linked to 
the finite-size induced features of the-J-V characteristics 
in the vicinity of the KT transitionrj 

Is the scaling property for the finite-size induced fea- 
ture of the I-V characteristics found in the present paper 
a real phenomena or just some accidental coincidence? 
Such accidental coincidences can of course never be com- 
pletely ruled out. However, our analysis really leaves 
little room for any such accidental coincidence. Yet, we 
admit that we do not have a plausible explanation at 
the moment. Thus the enigma with the scaling and the 
"ghost" transition calls for further theoretical, experi- 
mental and computational efforts. We have in the present 
paper described how the scaling property is reflected in 
the I-V data from a practical point of view and how to 
analyze the data from such a perspective. 

Finally we stress that our interpretation in terms of 
a finite-size induced feature of the I-V characteristics is 
entirely compatible with the established view of the KT 
transition with z — 2 and the usual finite-size scaling at 
a fixed T in the low-temperature phase, as we explicitly 
showed for the 2D RSJ model. In our finite-size inter- 
pretation there is no real phase transition but only a 
reflection of a "ghost" transition. This also agrees with 
the view by Strachan et al. in Ref. |3| that the apparent 
vanishing of the resistance in Recapi's data is not real but 
an artifact of the finite- voltage threshold in the measure- 
ments. On the other hand our conclusions are entirely 
incompatible with the interpretation put forward by Pier- 
son et at.Lj that there is a real resistive transition and 
that this is a KT transition with z « 6. Yet our analysis 
from a pragmatic poiat of view supports the original find- 
ing by Pierson et a^.llil that there is after all an intriguing 
scaling feature in the I-V characteristics associated with 
2D vortex fluctuations and that this scaling phenomena 
appears to be quite general. 
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Pierson et al. in Ref. [Ill argue, on the basis of a crude 
estimate of the finite size effects, that such effects cannot 
explain the Repaci data in contrary to the conclusion orig- 
inally reached by Repaci et al. in Ref. h2. In the present 
work we argue, on the basis of the strong similarity with the 
finite size induced scaling properties of the 2D RSJ model, 
that the Repaci data for small currents and lower temper- 
atures is finite size induced in agreement with Ref. n_2. 




FIG. 1. I-V characteristics for the 2D RSJ model on a 
square lattice with the finite size L = 8 plotted as _E = V/ L 
against J = I /L in log scales at temperatures T — 2.0, 1.5, 
1.3, 1.2, 1.1, 1.0, 0.95, 0.90, 0.85, 0.80, 0.75, 0.70, and 0.65 
(from top to bottom). The thick solid line corresponds to 
V oc /^ and has to a good approximation the same slope as 
the T — 0.90 curve over a finite-current region. For smaller 
currents the system becomes Ohmic, E = RJ, due to the 
finite size, which corresponds to the thick broken line. 
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FIG. 2. The I-V data for the 2D RSJ model with L = 8 
plotted as d\nE/dln J against J in log scale. From bottom to 
top the curves correspond to T = 2.0, 1.5, 1.3, 1.2, 1.1, 1.0, 
0.95, 0.90, 0.85, 0.80, 0.75, 0.70, and 0.65. The full drawn 
line segments connect the maxima. Down to d In E/d In J ~ 2 
the maxima occur at about the same J (between the vertical 
broken lines) and then move to higher J. The inset illustrates 
that the maxima of d In E/d In J to a good approximation are 
given by z-l- 1 where z is determined in two different ways (see 
text). The horizontal broken line in the inset corresponds to 
the KT value = 2 + 1 = 3 and the crossing point gives the 
estimate Tkt ~ 0.9. 
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FIG. 3. The size scaling given by Eq. (y) for the 2D RSJ 
model. L{E/ JY'^^'^' is plotted against LJ for a fixed T and 
the three different sizes L = 4, 6, and 8. The exponent z 
is determined from a scaling collapse of the data: (a) is for 
T = 0.90 and (b) for T = 0.80. In (a) and (b) a possible 
additional small L-dependence in Eq. (M) coming from ql is 
ignored and the values z ~ 2 and 3.3 are obtained for T = 0.90 
and 0.80, respectively. In (c) the L-dependence from ql is in- 
cluded for T = 0.80 and a good data collapse is obtained for 
the expected value z = 1/eT — 2 « 3.46 (see text). 
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FIG. 5. The function gL{T) determined for the 2D RSJ 
model with size L = 6 (open circles) and 8 (open squares). It 
is shown that the function ql for both sizes over a limited T 
interval is well approximated by ghiT) oc _R~" with a ~ 1/6. 
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FIG. 4. Demonstration of the existence of a scaling in the 
form of Eq. (M) for the size-dependent part of the I-V char- 
acteristics for the 2D RSJ model. E/JR is plotted against 
LJgLiT) and QhiT) is determined from the condition of a 
data collapse. In (a) QhiT) is determined for a fixed L for 
the two sizes L = 6 and 8 within the temperature intervals 
0.7 < T < 1 and 0.65 < T < 1, respectively. A good data 
collapse is obtained towards smaller values of LJgL{T) and 
furthermore the data for L = 6 and 8 also fall on top as re- 
quired by Eq. (W). In more detail, in (b) the same scaling 
plot is demonstrated for T — 0.9 (open squares), 0.8 (open 
circles), 0.7 (filled squares) and 0.65 (filled squares) in case of 
L = 8. 




FIG. 6. The function ghiT) as a function of z{T) for the 
2D RSJ model with size L = 8. The proportionality between 
gL(r) and R-'^^'^ holds in the interval 1.3 < z(T) < 3.5 (see 
Fig. I). 
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FIG. 8. The I-V characteristics for Repaci's YBCO data: 
V is plotted as a function of I for fixed T where 10K< T <40K 
(from left to right). The thick full line in the middle has the 
slope 3, which is the KT transition value and corresponds to 
the steepest slope of the data for T = 27K. The right thick 
full line has the slope 7, which corresponds to the "ghost" 
transition for the finite-size 2D RSJ model and which is the 
steepest slope of the data for T = 18K. The thick broken line 
has the slope 1 and corresponds to Ohmic resistance V = RI. 
Note the striking similarity with the 2D RSJ data in Fig. hi. 



FIG. 7. (a) Demonstration that the scaling obeyed by the 
finite-size induced part of the I-V characteristics for the 2D 
RSJ model with size L = 8 can be recast into the Pierson 
form: li+"/(r''°V") is plotted against IR-°'/T'^° and a 
good scaling collapse is found with a = 1/6 for the data 
in the interval 0.80 < T < 1.2. Identifying 1/a with z 
and R'" with ^ gives the connection R~^'^ ex ^ and makes 
the scaling look like the Pierson form. The inset shows 
that ln_R oc 1/VT'^'^ — 0.125 suggesting that the resistance 
would vanish at T^'^ = 0.125 and at this T (« 0.64) the 
dynamic critical exponent z(T) is to a good approximation 
given by z{T) — 1/T —2 = 6 consistent with the condi- 
tion z{T) = 1/a. (b) Pierson scaling plot including the lower 
temperatures T = 0.70 (open circles) and 0.65 (asterisks) for 
which ql ~ R"^'^' fails (see Fig. ]Sj. As seen in the upper 
left inset the same approximation for R can be used provided 

The lower right 



^yCG_ 0.125 is replac ed by VT^^ 
inset shows that \/T'^^ 



0.07. 



0.125 fails for the lowest values of T. 
In the main part, it is shown that a good scaling plot can still 
be obtained provided the two highest temperatures are ex- 
cluded but the ensuing z value is now z ~ 12. This illU|Hfcrates 
that the Pierson construction has too much flexibilityJlj 




FIG. 9. The I-V data for Repaci's YBCO data plotted 
as dlnV/dlnl against I in log scale for 17K< T <40K. The 
values of the maxima increase with decreasing T and give es- 
timates of z{T) -1-1. The horizontal broken line corresponds to 
the KT value z{T) + 1 = 3. The vertical broken line gives the 
approximate position for the maxima when dlnV/dlnl > 2, 
showing that the maxima occur at approximately the same 
/. The maxima for d\nV/dlnI < 2 occur at higher I with 
increasing T. Note the striking similarity with Fig. H for the 
finite-size 2D RSJ model. 
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FIG. 10. Demonstration of the existence of a scaling in the 
form of Eq. (|7[) for Repaci's I-V data. V/IR is plotted against 
Iql (T) and ql (T) is determined from the condition of a data 
collapse: (a) shows that a good scaling collapse is obtained 
for the data in the interval 23K< T <40K (the highest T 
corresponds to the lowest curve). Note the striking similarity 
with the 2D RSJ data in Fig. Wa). This similarity is further 
emphasized in (b) where the scaling function for the 2D RSJ 
model (full lines, note that E/J = V/I) is directly compared 
with the scaling function for the Repaci data (dotted lines). 
The two scaling function appears to fall on to each other when 
allowing for a constant shift along the x axis, suggesting that 
they are closely related. 



FIG. 12. Repaci's data analyzed according to Eq. (fel): 
<7l(T) in log scale is plotted against az{T)/[l — az{T)] where 
a according to Fig. hll is 1/6 to a good approximation and 
z{T) is determined from the maxima in Fig. y. As seen, the 
data fall approximately along a straight line from about where 
z{T) = 2 (vertical broken line). This corresponds to the T 
region where the maxima in Fig. occur for the same cur- 
rent I suggesting that this is the finite-size dominated part 
of the I-V characteristics. If the data would continue to fall 
on the straight line all the way to infinity, this would im- 
ply (/l(T') = oo and equivalently a vanishing resistance for 
z{T) = 1/q. Since this does not happen in practice we use 
the term "ghost" transition. 




FIG. 11. The function Pi(T') determined for Repaci's 
YBCO data: The filled circles give the values of ql (T) deter- 
mined from the data collapse in Fig. |lO|(a). The open circles 
demonstrates that the obtained QhiT) is well represented by 
3l(T) = AR-" with a = 1/6. Compare Fig. || for the 2D 
RSJ model. 




FIG. 13. Direct determination of the exponent a for 
Repaci's YBCO data. Eq. (X^ together with gi^ oc R~°' means 
that V/IR is a function of IR~°' . V{= RI) is plotted as a 
function of / for a given fixed value of V/IR =const. (open 
and filled circles correspond to const. =1.5 and 2.0, respec- 
tively). Consequently, the scaling implies that the data points 
should fall on straight lines given by RI oc j-^i°'-^ and the 
slope of these lines should be given by —1/a — 1. As seen this 
prediction is borne out and the obtained values are a~^ = 6.3 
and 6.2 in close agreement with the value a~^ ~ 6 obtained 
from Fig. 0. 
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Fip^ 14. The Pierson scaling plot for Repaci's YBCO 
datallll 7i+i/7(rFi/^) is plotted against I^/T where it is 
assumed that ^ oc R~ '^ . The data is plotted in log scales 
and the construction makes the scales very extended which 
exaggerates the goodness of the data collapse. The high-T 
branch of the scaUng plot corresponds to 18K< T <40K and 
is given by the lower scaling curve and should be compared 
to the scaling of the finite-size 2D RSJ model in Fig. W. The 
low-T branch contains data for 10K< T <16K. A transition 
of a vanishing resistance at T = 17K jwith z — 5.9 was con- 
cluded from this scaling construction.tJ We claim that this 
scaling plot does not reflect any real transition to a vanishing 
resistance. 
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